Bromophenyl functionalization of carbon nanotubes : an ab initio study. 
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We study the thermodynamics of bromophenyl functionalization of carbon nanotubes with respect 
to diameter and metallic/insulating character using density-functional theory (DFT). On one hand, 
we show that the activation energy for the grafting of a bromophenyl molecule onto a semiconducting 
zigzag nanotube ranges from 0.73 eV to 0.76 eV without any clear trend with respect to diameter 
within numerical accuracy. On the other hand, the binding energy of a single bromophenyl molecule 
shows a clear diameter dependence and ranges from 1.51 eV for a (8,0) zigzag nanotube to 0.83 eV 
for a (20,0) zigzag nanotube. This is in part explained by the transition from sp^ to sp'^ bonding 
occurring to a carbon atom of a nanotube when a phenyl is grafted to it and the fact that smaller 
nanotubes are closer to a sp^ hybridization than larger ones due to increased curvature. Since a 
second bromophenyl unit can attach without energy barrier next to an isolated grafted unit, they 
are assumed to exist in pairs. The para configuration is found to be favored for the pairs and their 
binding energy decreases with increasing diameter, ranging from 4.34 eV for a (7,0) nanotube to 
2.27 eV for a (29,0) nanotube. An analytic form for this radius dependence is derived using a tight 
binding hamiltonian and first order perturbation theory. The dependance obtained (where R 

is the nanotube radius) is verified by our DFT results within numerical accuracy. Finally, metallic 
nanotubes are found to be more reactive than semiconducting nanotubes, a feature that can be 
explained by a non-zero density of states at the Fermi level for metallic nanotubes. 



I. INTRODUCTION 

Since their discoverjW^I, a lot of efforts have been de- 
voted to harnessing the exceptional characteristics of 
single-walled carbon nanotubes into useful electronic de- 
vices. Despite numerous realizations of such devices,'^lEl 
schemes for their production at the commercial scale are 
scarce. The main problem lies in the selection and ma- 
nipulation of a single nanotubes, mainly because the van 
der Waals force causes them to agglomerate in bundles 
Moreover, the high molecular weight and hydrophobic 
character of carbon nanotubes makes them particularly 
difficult to s olubilize and to handle on a substrate or in 
solution.laESl 

A proposed solution is to functionalize covalently the 
nanotube walls to modify their interaction with their 
environment.'^ A high yield method developed by J.M. 
Tour et al. uses the grafting of a phenyl unit.C^ These 
are commonly added to different molecules to ease their 
manipulatiorP^ and increase their solubility.'^ It is fur- 
thermore reversible,'^ allowing the original electronic 
features of the nanotubes to be retrieved by thermal 
annealing.'i^ Although the mechanism governing this re- 
action has been previously described in the literature, 1- 
the dependence of its thermodynamics on the nanotube 
diameter and metallicity is not fully known and is there- 
fore investigated in this work. 

To be more specific, we report ab initio results on 
the thermodynamics of the functionalization of carbon 
nanotubes with bromophenyl molecules. We study bro- 
mophenyl since it has more practical interest that phenyl, 
that is, because the bro mine atom can be used in a 
successive Suzuki reactiorP^E^ to attach a vast number 



of other organic molecules to the grafted phenyl. This 
scheme is commonly used, for instance, in nanotube het- 
eroj unction design.l^l 

This article is organized as follows : the next section 
describes the computational procedure used, followed by 
results on the thermodynamics of the functionalization 
of a carbon nanotube by a single bromophenyl unit. 
Next, we assess the thermodynamics of pair formation 
by grafted bromophenyls. Finally, we derive an analyti- 
cal form for the diameter dependance of the binding en- 
ergy from a tight binding hamiltonian and first order 
perturbation theory, then verify its agreement with our 
computed results. 



II. COMPUTATIONAL DETAILS 

All the results in this work wer e produced using 
density-functional theory (DFT)?^^ as implemented in 
the ONETEF^^ code. The reduced computational ef- 
fort for this implementation of DFT (linear scaling when 
applied to semiconductors) enables us to investigate the 
bromophenyl functionalization for a wider range of sys- 
tem sizes. Also, the use of Nonorthogonal Generalized 
Wannier Functions (NGWFspSJ ensures that the amount 
of empty space in the simulated system has little im- 
pact on the computation time, thus easing the simula- 
tion of large nanotubes. The system simulated in the 
present work consists of five zigzag nanotube primitive 
cells, with a total length of 2.13 nm and 0, 1 or 2 bro- 
mophenyl units functionalized on it. The chirality of the 
nanotubes is chosen to be zigzag since it keeps the primi- 
tive cell size minimal and allows us to simulate both semi- 
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conducting and metallic nanotubes. Periodic boundary 
conditions are applied to this system, so that the effec- 
tive length of the simulated nanotube is infinite. Putting 
one single bromophenyl or bromophenyl pair per 5 zigzag 
nanotube primitive cells ensures that periodic images of 
the bromophenyl functionalization will not interact with 
each other. Therefore, this system size provides fully 
converged binding and activation energies, even for the 
metallic nanotubes. Due to the reduced computational 
effort required, we were able to investigate nanotubes of 
chirality up to (29,0). Also, the nanotubes are kept 2 nm 
apart from their periodic images in the direction perpen- 
dicular to the nanotube axis, to prevent any spurious 
nanotube-nanotube interaction. Since the nanotubes are 
uncharged and have no large dipole moment, truncation 
of the Coulomb interactiorP^ is not required in this case, 
but could be used for other ligand species. 

The electronic structure is represented by 1, 4 and 7 
NGWFs with a cutoff radius of 8.0, 7.0 and 8.0 bohr 
for the hydrogen, carbon and bromine atoms respec- 
tively. An effective cutoff energy, analogous to the one in 
plane wave pseudopotential codes, is set to 35 Ha, which 
is quite sufficient for our systems (c.f. E.R. Margine 
and . The PBEpS approximation for the exchange- 
correlation energy is used throughout this work. Spin- 
polarized calculations are used whenever the system con- 
tains an odd number of electrons. A level of convergence 
was achieved such that the calculated activation and 
binding energies have a numerical accuracy of 50 meV 
per primitive cell (consisting of 5 zigzag nanotube prim- 
itive cells). Geometric relaxations of each structure were 
then carried out with a tighter convergence criteria of 
10 meV per primitive cell. The optimization of the 
NGWFs within ONETEP results in a smooth poten- 
tial energy surface from which accurate ionic forces can 
be calculatecP^ and eliminates the basis set superposi- 
tion erroi^ present in traditional DFT implementations 
based on localized basis sets. This feature allows conver- 
gence criteria as tight as the one above to be reached. 



III. ACTIVATION ENERGY 
A. Procedure 

In this section, we report our calculated activation and 
binding energies for the functionalization of a nanotube 
with one bromophenyl unit. In both cases, the first step 
is to relax a system where the pristine nanot ube and 
the bromophenyl unit are 7 A apart (see Fig. 1(a)) to 
find the ground state energy of the isolated reactants. 
Then, a functionalized nanotube (see Fig. 1(b) I is relaxed 
to obtain the ground state energy of the product. The 
binding energy is then obtained by taking the difference 
between the two results. 

To find the activation energy, we linearly interpolate 
the atomic positions of 12 intermediate systems between 
the isolated reactants and the final product. The ground 




(a) Initial configuration (b) Final configuration 

FIG. 1. Initial and final configurations used to simulate the 
bromophenyl functionalization of a (9,0) carbon nanotube. 
The distance between the nanotube and the bromophenyl unit 
is set to 7 A in the initial configuration to prevent any inter- 
action between the two systems. In the final configuration, 
the relaxed bromophenyl-nanotube bond length is 1.54 A. 



state energy of each interpolated system is then calcu- 
lated without any geometric relaxation. The top of the 
energy barrier is taken to be the highest energy value 
obtained among the 12 interpolated calculations. The 
activation energy is obtained by taking the difference 
between the top of the energy barrier and the ground 
state energy of the isolated reactants. As an example, a 
graph of the total energy with respect to bromophenyl- 
nanotube distance is given for a (13,0) nanotube in Fig. [2] 
We obtain an activation energy of 0.76 eV and a binding 
energy of 1.06 eV which are of the right order of magni- 
tude according to previous studies.l^ 

It should be noted that the activation energies ob- 
tained here are in fact upper bounds to the exact DFT 
ones. Rigorously, a complete relaxation for each geomet- 
ric distance between the nanotube and the bromophenyl 
unit should have been performed. This would have low- 
ered slightly our estimation of the activation energy and 
greatly increased the computational cost. Still, this 
should not prevent diameter or metallicity dependence 
from being observed in our calculations. However, our ac- 
tivation energy calculations might overestimate by more 
than 50 meV the exact DFT prediction, despite the tight 
numerical precision of our calculations. 

To get an intuitive understanding of the activation and 
binding energies trends with respect to nanotube diam- 
eter and metallicity, we divided both kind of energies 
into three parts : the contribution from the deformation 
of the nanotube, the contribution from the deformation 
of the bromophenyl unit and the remainder, labeled the 
electronic contribution. For the binding energy, the nan- 
otube deformation contribution is obtained from the dif- 
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B. Results 
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FIG. 2. Reaction barrier for the functionalization of the 
(13,0) nanotube by a bromophenyl unit. The dashed line is 
a guide for the eye and is set to the ground state energy of 
the reactants. Each point between the initial and final ones 
is a linear interpolation between the two geometries shown in 



TABLE I. Contributions to activation energies (in eV) for a 
bromophenyl unit functionalization. 





Semiconducting 


Metallic 


Nanotube 


(8,0) 


(13,0) 


(20,0) 


(9,0) 


Nanotube deformation 


0.70 


0.81 


0.86 


0.63 


Bromophenyl deformation 


0.06 


0.06 


0.06 


0.10 


Electronic contribution 


-0.03 


-0.11 


-0.15 


-0.09 


Activation energy 


0.73 


0.76 


0.77 


0.64 



ference between the total energy of a system having the 
same geometry as the functionalized nanotube, but from 
which the bromophenyl has been removed, and the total 
energy of the pristine nanotube. The same scheme was 
applied for the activation energy, except that the func- 
tionalized nanotube is replaced with the system configu- 
ration at the top of the energy barrier. The bromophenyl 
deformation contribution is obtained with a scheme anal- 
ogous to the nanotube deformation contribution. 



TABLE II. Contributions to binding energies (in eV) for a 
bromophenyl unit functionalization. 





Semiconducting 


Metallic 


Nanotube 


(8,0) 


(13,0) 


(20,0) 


(9,0) 


Nanotube deformation 


1.12 


1.30 


1.38 


1.17 


Bromophenyl deformation 


0.11 


0.12 


0.11 


0.12 


Electronic contribution 


-2.74 


-2.47 


-2.33 


-3.14 


Binding energy (absolute value) 


1.51 


1.05 


0.84 


1.85 



Results for the activation energy are presented in Ta- 
ble|lj For the semiconducting nanotubes, the deformation 
contribution increases with their diameter. However, the 
electronic contribution has the opposite trend, so that the 
two effects cancel out up to numerical precision. Thus, 
the activation energy is found to be between 0.73 and 
0.76 eV, without any significant diameter dependence. 
However, the metallic nanotube (9,0), with an activation 
energy of 0.64 eV, shows a reactivity significantly higher 
than the semiconducting nanotubes. This trend agrees 
with other DFT studies on NOj^SI and carben^^U nan- 
otube functionalization and is associated with a non-zero 
electronic density of states (DOS) at the Fermi level of 
metallic nanotubes.!^ As we shall see later, all metallic 
nanotubes in the present study are more reactive than 
semiconducting tubes of similar diameter. 

Results for the binding energy are presented in Ta- 
bic [nj For the semiconducting nanotubes, the deforma- 
tion energy cost increases with diameter. This is easy 
to understand intuitively. The bromophenyl functional- 
ization involves changing the hybridization of a carbon 
atom of the sidewall of the nanotube from sp^-like to 
sp^-M^e.^'^ For larger diameters, the geometry of the wall 
is close to planarity and to ideal geometry of sp^ bond- 
ing. Thus, when the nanotube's diameter increases, it 
becomes more energetically unfavorable to undergo geo- 
metrical deformation to sp^ bonding. That explains why 
the deformation contribution to the binding energy fa- 
vors functionalization of small nanotubes over large ones. 
This analysis applies both to binding and activation ener- 
gies, thus explaining why deformation energy favors small 
tube functionalization in both cases. 

The electronic contribution to the binding energy also 
favors small nanotubes. Since the bromophenyl deforma- 
tion energy shows no trend with respect to diameter and 
both other contributions favor small nanotubes, the total 
binding energy is stronger when the nanotube is smaller. 

The binding energy trend with respect to metallicity 
shows, for the same reason given in the activation en- 
ergy discussion, that the electronic contribution is higher 
for the metallic nanotubes than for the semiconducting 
ones. Since both deformation contributions to the bind- 
ing energy remain similar for nanotubes of similar size, 
whether they are semiconducting or metallic, the over- 
all result is that metallic nanotubes bond more strongly 
with bromophenyl than semiconducting ones. 

In summary, our results show that the activation en- 
ergy of bromophenyl functionalization of semiconducting 
nanotubes has almost no diameter dependence, although 
the binding energy is significantly higher for small tubes, 
making their functionalization more stable. This is in 
part explained by the fact that nanotubes with higher 
curvature are closer to a sp'^ hybridization. Also, thanks 
to the presence of a non-zero DOS at their Fermi level, 
metallic nanotubes can be functionalized by bromophenyl 
units with more ease and more stability than semicon- 
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FIG. 3. Energy levels and their corresponding electronic den- 
sities for orbitals near the Fermi level of functionalized (13,0) 
nanotube. Only the functionalized site and the surrounding 
atoms are shown. 

ducting nanotubes of similar diameter. This could pro- 
vide an easy way to sort semiconducting nanotubes from 
metallic ones. 



IV. FUNCTIONALIZATION IN PAIRS 
A. Procedure 

The study presented in the previous section explores 
how a bromophenyl unit reacts with a pristine nanotube. 
However, it has been shown in the literature that, at 
room temperature, a single bromophenyl grafted onto a 
nanotube can desorb or migrate on its surface.fi^ There- 
fore, a stable bromophenyl unit on a nanotube cannot be 
isolated as pictured in Fig. [l] It is believed that it must 
instead be adjacent to another bromophenyl unit.^^ The 
motivation for this hypothesis can be explained qualita- 
tively and was introduced by the groups of C. Dyke'^-'^ and 
G. Schmidt-ISS] Just before grafting onto the nanotube, the 
bromophenyl radical gets an unpaired electron from the 
separation of phenyldiazonium (i?C6H4Nj) and is there- 
fore highly reactive. 

When the first bromophenyl unit covalently bonds, its 
half-filled orbital mixes with the orbitals of the tube. 
This adds two new electronic states in the gap of the 
pristine nanotube according to our calculations. One 
is associated with the filled state of the half-filled or- 
bital of the bromophenyl unit and constitutes the high- 
est occupied molecular orbital (HOMO) of the function- 
alized nanotube. The other electronic state is associ- 
ated with the empty state of the half-filled orbital of 
the bromophenyl unit and constitutes the lowest unoc- 
cupied molecular orbital (LUMO) of the functionalized 
tube. Fig. |3] shows the electronic density of the orbitals 



near the Fermi level of the functionalized (13,0) nan- 
otube. Only a section of the sidewall of the tube close 
to the functionalized site is shown to ease visualization. 
We notice that the LUMO has exactly the same shape as 
the HOMO state and is 50 meV above the Fermi level. 
Both the HOMO and LUMO states are weh localized 
near the functionalized site, whereas the other orbitals 
are spread over the nanotube. The other nearest states 
are further apart in energy from the Fermi level, about 
350 meV from it. Since the new LUMO only exists a few 
angstroms around the first bromophenyl unit on the nan- 
otube's sidewall, the functionalization of the second bro- 
mophenyl, eased by the higher electronic affinity around 
the first one, takes place close to it. Thus, at room tem- 
perature, bromophenyls should be found in pairs on the 
surface of a functionalized nanotube. 

Since we are interested in the thermodynamics of the 
functionalization, it is relevant to investigate the stability 
of the final product obtained, that is, the binding energy 
of pairs of bromophenyls on carbon nanotubes. Yet, a 
calculation of the binding energy of the pairs requires 
a precise knowledge of the configuration of these pairs. 
Since the LUMO is strongly localized around the first 
bromophenyl unit, we expect the second bromophenyl to 
bind more strongly closer to the first one. We then have 
three likely configurations, pictured in Fig. [4] and named 
ort/io(first adjacent atom), meto(second adjacent atom) 
and para (third adjacent atom). To evaluate the likeliness 
of the different reaction sites, let's first neglect the effect 
of the sidewall curvature of the nanotube. 

It has been shown in the literature, using ab initio 
calculations, that the meta configuration of aryl groups 
on a graphene sheet is unstable'^'^ in accordance with the 
aryl addition rule. Moreover, the LUMO state illustrated 
in Fig. [3] is localized on the ortho and para sites, but not 
the meta sites. The functionalization of the former ones 
is therefore favored over the latter. 

However, steric constraints will lower the binding en- 
ergy of the second phenyl if it binds on atoms adjacent 
to the first reaction site, thus disfavoring the ortho site. 

In addition, the binding energies of phenyls on 
graphene for these configurations were found by Margine 
et al. to be 1.25 eV for the ortho site and 1.51 eV for the 
para site.'^ The higher binding energy of the para config- 
uration over the ortho configuration shows that although 
the ortho configuration has larger LUMO electronic den- 
sity than the para configuration, the steric constraints 
disfavor the closest site. 

In our study, a hydrogen of the phenyl is replaced by 
a bromine atom, which should further increase the steric 
constraints. Indeed, our results also show that para con- 
figuration is favored over the ortho configuration in our 
systems. Therefore, we limit our study to the former 
case. 

The reader should note that no reaction barriers are 
shown for the second functionalization. First, it is very 
difficult to define the reaction pathway at this stage. Sec- 
ond, the activation energy for the second bromophenyl is 
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FIG. 4. Three possible configurations for a bromophenyl 
pair. The "x" shows the first functionalized site. The location 
of the second unit gives the ortho (1), meta (2) and para (3) 
configurations. 



known to be smaller than for the first reaction^^ and thus 
does not limit the reaction rate. 

If we now take into account the curvature of the side- 
wall, we find that there are two possible variants for the 
para configuration on a zigzag nanotube. The imaginary 
line linking the two grafted sites can either be parallel to 
the axis of the tube or make an angle of 60° with it. Both 
cases should show different binding energies since the cur- 
vature of the nanotube affects their geometry. Therefore, 
for pairs of bromophenyls in the para configuration, we 
investigate the magnitude of this difference in addition to 
the diameter and metallicity dependence of their binding 
energy. 
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FIG. 5. Binding energy of bromophenyl pairs on zigzag 
nanotubes for two para variants (0° and 60° relative to 
the tube axis). The continuous line shows the best fit of 
Eq. |6] on the semiconducting nanotubes data. We obtain 
-Bbinding = 15.00 eV k^/R^ + 2.11 eV. The standard devi- 
ation a = 86 meV is similar to our numerical accuracy of 
50 meV, which supports the validity of our model (Eq. [6|. 



TABLE III. Binding energies (in eV) for bromophenyl pairs 
{para configuration, 0° relative to the axis, see Fig.[5| and for 
two isolated units. 





Semiconducting 


Metalhc 


Nanotube 


(8,0) 


(13,0) 


(20,0) 


(9,0) 


Bromophenyl pair 


3.63 


2.76 


2.38 


3.77 


Single bromophenyl unit x 2 


3.02 


2.12 


1.66 


3.72 



B. Results 

The results are shown in Fig. [sj The difference of bind- 
ing energy between the two variants of para configura- 
tions (0° and 60°) does not show any clear trend with 
respect to diameter or metallicity and its magnitude re- 
mains small in all cases. Therefore, we cannot conclude 
whether there is a variant favored in the functionalization 
process. 

However, the binding energy trend with respect to 
diameter is clear : bromophenyl pairs are more stable 
on small nanotubes. This conclusion agrees with previ- 
ous ab initio studies of nanotube functionalization with 
phenyls^^ and N02.'^^ The fit in Fig.[5]is explained in the 
next section. 

Table |III| compares the binding energy of para paired 
bromophenyls units and isolated units. The results con- 
firm the pair configuration to be more stable than the 
isolated configuration for all semiconducting cases stud- 
ied. However, they cast doubt in the case of metallic 
nanotubes since the pairs are only 50 meV more stable 
than isolated units and could therefore separate easily at 



room temperature (fc^T = 25 meV). This observation is 
in agreement with the explanation given earlier for the 
pairing of bromophenyls. Indeed, since the pairing occurs 
because the LUMO of a nanotube functionalized with a 
single bromophenyl is well inside the electronic gap of 
the pristine nanotube, the proposed mechanism cannot 
account for an eventual pairing of the bromophenyls on 
metallic nanotubes since they have no gap. 



V. THEORETICAL BACKGROUND 

It is interesting to derive an approximate expression 
for the diameter dependance of the binding energy of 
a bromophenyl pair on a nanotube. At first order in 
perturbation theory, it is given by the expectation value 
of the Hamiltonian of the bromophenyl pair, iJ', in the 
eigenstate of the pristine nanotube (most) involved in the 
bond formation, | 7r„). Since the other eigenstates have 
a smaller overlap with the orbitals of the phenyl graft, 
we do not expect them to contribute significantly to the 
bromophenyl-nanotube bond. 
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With a tight binding approach, following the notation 
of M. Pudlak et a/.f^the |7r) orbital of a nanotube can 
be written as: 



, , , /3cos20 - 1 , , 

K„)=tan0^/ 3^^^^^ |.) 



tan e . , Vcos 29 , , 



(1) 



%/3 cos 26 ' ' cos 
where is a measure of the curvature and is defined as: 

• " (2) 



where a = 1.42 A is the length of a covalent bond in the 
nanotube and R is the radius. To first order in static 
non-degenerate perturbation theory, we have: 



{s\H'\py) = {s\H'\p,) = {py\H'\p,) = 0. (3) 



Therefore, 



> 9 ^/3cos2(? — 1\ , . , 

^„ i?>„ ^tan^e — — {s H' s) 

V 3 cos 20 / 



tan^ 



3 cos 261 
cos 29 



{Pv\H'\Py) 



(4) 



cos^ 



{Pz\H'\p.) 



With a Taylor expansion and some algebra, we find that 
for small values of 1/R (i? ^ a), at leading order : 



with a standard deviation of cr = 86 meV. However, for 
metallic nanotubes, the amount of available data was in- 
sufhcient for a fit of Eq. [6] to be meaningful. 

The discrepancy between the fit on semiconducting 
nanotube data (Eq. uh and our results is similar to our 
numerical accuracy [crtit — 86 meV ^ cTconvorgcnco ~ 
50 meV). Also, our fitted value of c is 2.11 eV, in 
very close agreement to the calculated binding energy 
of 2.10 eV, which is not part of the fitted data. There- 
fore, Eq. [6] accurately describes the diameter dependance 
of semiconducting nanotubes functionalized with bro- 
mophenyl pairs. Since no properties exclusive to bro- 
mophenyl has been assumed in the derivation, it should 
accurately describe any covalent functionalization of nan- 
otubes where first order perturbation theory remains 
valid, i.e. where the binding energy is of the same order 
of magnitude as the cases studied here. 

The binding energy also shows a clear trend with 
metallicity : it is significantly higher for metallic nan- 
otubes than for semiconducting nanotubes of similar di- 
ameter. This is expected since the explanation given for 
the activation energy trend with metallicity remains valid 
for the binding energy of pairs. Together with the acti- 
vation energies, these results confirm that bromophenyl 
functionalization thermodynamically favors nanotubes of 
smaller diameter and metallic nanotubes over semicon- 
ducting ones. 



VI. CONCLUSION 



(^„|i?>„) = (p^H'Ip,) + ^- {s\H'\s) 

+ l{Py\H'\Py)~{p.\H'\p^)']^. 



(5) 



Therefore, the trend with respect to diameter of the 
binding energy for a bromophenyl pair functionalization 
can be written as: 



-^binding — ^ ^2 ~^ 



(6) 



where c is the binding energy for graphene. This argu- 
ment also applies to the binding energy of a single bro- 
mophenyl. However, in the latter case, we do not have 
enough results to assess the accuracy of Eq. [6j 

A fit of Eq. [6] on binding energies for pairs of bro- 
mophenyl on semiconducting nanotubes is shown in 
Fig. [5] The agreement is excellent. We find: 

A' 

i^binding = 15.00 eV-^ + 2.11 eV, (7) 



Our activation and binding energy calculations show 
that bromophenyl functionalization is selective with re- 
spect to nanotube size and metallicity : functionalization 
of smaller and/or metallic nanotubes is thermodynami- 
cally favored. However, this selectivity is stronger for 
smaller nanotubes while the experimental tubes are rela- 
tively large. Therefore, it may be challenging to observe 
this selectivity in experimental samples, due to others 
factors becoming of equal importance as compared to the 
diameter, such as the family of the nanotube. 
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